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1. Introduction 

Recently, spatially homogeneous imperfect fluid models have been investigated using 
techniques from dynamical systems theory jl], [^, |3|, ffl. In these papers, dimensionless 
variables and a set of dimensionless equations of state were employed to analyze various 
spatially homogeneous imperfect fluid cosmological models. It was also assumed that 
the fluid is moving orthogonal to the homogeneous spatial hyper surf aces; that is, the 
fluid 4-velocity, u a , is equal to the unit normal of the spatial hypersurfaces. The energy- 
momentum tensor can be decomposed with respect to u a according to f|: 

T ab = (p + P)u a u b + pg ab + q a u b + u a q b + n ab , (1.1) 

where p = p + II and p is the energy density, p is the thermodynamic pressure, II is the 
bulk viscous pressure, 7r ab is the anisotropic stress, and q a is the heat conduction vector 
as measured by an observer moving with the fluid. 

In papers [|, |J it was also assumed that there is a linear relationship between 
the bulk viscous pressure II and the expansion 9 (of the model), and a linear relationship 
between the anisotropic stress 7r a b and the shear a ab ; that is, 

IT = - CO (1.2a) 
n ab = - 2r]a ab (1.26) 

where ( denotes the bulk viscosity coefficient and 77 denotes the shear viscosity 
coefficient. Equations (1.2a) and (1-26) describe Eckart's theory of irreversible 
thermodynamics. Eckart's theory is a first order approximation of the viscous pressure 
II and the anisotropic stress n ab and is assumed to be valid near equilibrium 0. Coley 
and van den Hoogen ]IJ and Coley and Dunn || have studied the Bianchi type V 
cosmological models using equations (1.2) as an approximation for the bulk viscous 
pressure and the anisotropic stress. They found that if the models satisfied the weak 
energy conditions, then the models necessarily isotropize to the future. Belinskii and 
Khalatnikov |7j, with different assumptions on the equations of state, found similar 
behaviour present in the Bianchi type I models. The addition of viscosity allowed for 
a variety of different qualitative behaviours (different from that of the corresponding 
perfect fluid models). However, since the models studied in |], 0, ||, [| satisfy Eckart's 
theory of irreversible thermodynamics ||, they suffer from the property that signals in 
the fluid can propagate faster than the speed of light (i.e., non-causality), and also that 
the equilibrium states in this theory are unstable (see Hiscock and Salmonson || and 
references therein). Therefore, a more complete theory of irreversible thermodynamics 
is necessary for fully analyzing cosmological models with viscosity. 

Among the first to study irreversible thermodynamics were Israel and Israel and 
Stewart [[10], |TT|. They included additional linear terms in the relational equations (1.2). 



Assuming that the universe can be modelled as a simple fluid and omitting certain 
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divergence terms, the linear relational equations for the bulk viscous pressure, the heat 



conduction vector, and the shear viscous stress are [J I 



n = -C«« + A)n-a g? B ), (1.3a) 
q - = - K Th a \T~ 1 T. b + u b + p 1 q b -a U. b -a 1 n% c ), (1.3b) 
7t ab = -2r](u a;b + (3 2 Tx a b- onq a;b ), (1.3c) 

where h ab = u a u b + g ab is the projection tensor and (A ab ) = \h c a h d b (A cd + A dc - 
\h C( ih e f A e f). The variable T represents the temperature, k represents the thermal 
conductivity, Po, Pi and Pi are proportional to the relaxation times, «o is a coupling 
parameter between the heat conduction and the bulk viscous pressure and «i is a 
coupling parameter between the shear viscous stress and the heat conduction. We shall 
refer to equations (1.3) as the truncated Israel-Stewart equations. These equations, 
(1.3), reduce to the Eckart equations (1.2) used in Jj], |^, |3|, £| when «o = ot\ = Po = 

Pi =lh = 0. 

Belinskii et al. [12|| were the first to study cosmological models satisfying the 



truncated Israel- Stewart theory of irreversible thermodynamics. Using qualitative 
analysis, Bianchi I models were investigated with a relational equation for the bulk 
viscous pressure and the shear viscous stress of the form (1.3). They also assumed 



equations of state of the form |L2 



C = Cop m , V = Vop n , Po = p-\ and P 2 = p~\ (1.4) 

where m and n are constants and ( and r] are parameters. The isotropizing effect 
found in the Eckart models no longer necessarily occurred in the truncated Israel- 
Stewart models. It was also found that the cosmological singularity still exists but 
is of a new type, namely one with an accumulated "visco-elastic" energy [|l!|. Similar 
to the work done by Belinskii et al. |~2|, Pavon et al. fl3j and Chimento and Jakubi 



Tj| studied the flat Friedmann-Robertson- Walker (FRW) models. They assumed the 



same equations of state as Belinskii et al. []12 , namely ( |1.4|) , but studied the models 



using slightly different techniques. Chimento and Jakubi |L4| also found exact solutions 
in the exceptional case m = 1/2 (which will be of interest later). They found that the 
future qualitative behaviour of the model was independent of the value of m; however, 
to the past, "bouncing solutions" and deflationary evolutions are possible [fl4" |. 



In addition, Hiscock and Salmonson || investigated further generalizations of the 
truncated Israel- Stewart theory of irreversible thermodynamics. Namely, they included 
the non-linear divergence terms in the relational equations (1.3). We shall refer to such 
a theory as the "Israel-Stewart-Hiscock" theory. Hiscock and Salmonson used equations 
of state arising from the assumption that the fluid could be modelled as a Boltzmann 
gas. They concluded that when the Eckart equations, (1.2), or the truncated Israel- 
Stewart equations, (1.3), were used, inflation could occur, but if the non-linear terms 
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of the Israel-Stewart-Hiscock theory were included, inflation was no longer present. 
This result led them to conclude that "inflation is a spurious effect produced from 
using a truncated theory" ||. However, Zakari and Jou |D| also employed equations 
arising from the full Israel-Stewart-Hiscock theory of irreversible thermodynamics, but 
assumed equations of state of the form (|1.4j) and found that inflation was present in all 
three theories (Eckart, truncated Israel-Stewart, Israel-Stewart-Hiscock). Therefore, it 
appears that the equations of state chosen determine if the model will experience bulk- 
viscous inflation. Romano and Pavon JTJJ also analyzed Bianchi III models using both 
the truncated Israel-Stewart theory and the Israel-Stewart-Hiscock theory. They only 
analyzed the isotropic singular points, but concluded that the qualitative behaviour of 
the models in the two different theories was similar in that the anisotropy of the models 
dies away and the de Sitter model is a stable attractor. 

In this work we use equations of state that are dimensionless, and hence we are 
generalizing the work in Jl|, ^ |3|, |J in which viscous fluid cosmological models satisfying 
equations (1.2) were studied. One reason for using dimensionless equations of state is 
that the equilibrium points of the system of differential equations describing spatially 
homogeneous models will represent self-similar cosmological models [17]]. In addition, it 
could be argued that the use of dimensionless equations of state is natural in the sense 
that the corresponding physics is scale invariant (see also arguments in Coley ||18|| ). 

The intent of this work is to build upon the foundation laid by Belinskii et 
al. [0, Pavon et al. |13| and Chimento and Jakubi and investigate viscous 



fluid cosmological models satisfying the linear relational equations (1.3). We will use 
dimensionless variables and dimensionless equations of state to study the qualitative 
properties of isotropic and spatially homogeneous cosmological models. In particular, 
we shall study this new "visco-elastic" singularity and we shall determine whether bulk- 
viscous inflation is possible. We will also determine if there is a qualitative difference 
between these models and the models studied by Burd and Coley (3| where the Eckart 
equation (1.2 a) was assumed. 

In section ^| we define the models and establish the resulting dynamical system. In 
section |3| we investigate the qualitative behaviour of the system for different values of the 
physical parameters. In section |4] we discuss and interpret our results and in section ^| 
we end with our conclusions. For simplicity we have chosen units in which 8ttG — c — 1. 



2. Friedmann- Robertson- Walker Models 

In this paper we assume that the spacetime is spatially homogeneous and isotropic and 
that the fluid is moving orthogonal to the spatial hypersurfaces. The energy-momentum 
tensor considered in this work is an imperfect fluid with a non-zero bulk viscosity (that 
is there is no heat conduction, q a = 0, and no anisotropic stress, ir ab = 0). The Einstein 
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field equations (G a b = T ab ) and the energy conservation equations {T ab b u a = 0) can be 
written as (see Burd and Coley |4[]): 

9 = - V-I( p + 3p), (2.1a) 

p = -9(p + p), (2.16) 

9 2 = 3p-^% (2.1c) 

where 9 is the expansion and 3 R is the curvature of the spatial hypersurfaces. If the 
curvature is negative, i.e., 3 R < 0, then the FRW model is open, if 3 R = then the 
model is flat, and if 3 R > then the FRW model is closed. Assuming that the energy 
density, p, is no n- negative, it is easily seen from (2.1c) that in the open and flat FRW 
models the expansion is always non-negative, i.e. 9 > 0, but for the closed FRW models 
the expansion may become negative. (Great care must be taken in this case because 
the dimensionless quantities that we will be using become ill-defined at 9 = 0.) 
We can obtain an evolution equation for II by solving (1.3a) for IT, 

* = (2 ' 2) 

Now the system of equations defined by equations (2.1a), (2.16), and Q2.2| ) constitute 
a dynamical system of the form X = F(X), where X = (#, p, IT). This system of 
equations is invariant under the mapping (see Coley and van den Hoogen [I], [T7|j ) 

0^x9, n->A 2 n, P ^x 2 p, 

P - A 2 P , C - AC, A, - a- 2 a,, t - \-H, 1 ' j 

and this invariance implies that there exists a symmetry in the dynamical system ||19|| . 
Therefore, we introduce new dimensionless variables x, y and a new time variable Q as 
follows: 

and ^ = -3^' (2 ' 4) 
and consequently the Raychaudhuri equation, (2.1a), effectively decouples from the 
system. 

In order to complete the system of equations we need to specify equations of state 

for the quantities p, Po, and (. In principle equations of state can be derived from 

kinetic theory, but in practice one must specify phenomenological equations of state 

which may or may not have any physical foundations. Following Coley |18|, ^0(1 , we 

introduce dimensionless equations of state of the form 
p_ 

ff 1 - 

~ = CoX m } (2.56) 
ax ri , (2.5c) 



= p x\ (2.5a) 
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where p Q , ( , and a are positive constants, and £, m, and r± are constant parameters (x is 
the dimensionless density parameter defined earlier). In the models under consideration, 
9 is positive in the open and flat FRW models, thus equations (2.6) are well defined. In 
the closed FRW model the expansion could become zero, in which case these equations 
of state break down. However, we can utilize these equations to model the asymptotic 
behaviour at early times, i.e., when 9 > 0. The most commonly used equation of state 
for the pressure is the barotropic equation of state p — (7 — l)p, whence p D — |(7 — 1) 
and I = 1 (where 1 < 7 < 2 is necessary for local mechanical stability and for the speed 
of sound in the fluid to be no greater than the speed of light). 

We define a new constant b = a/( Q . Using equations (2.5) and (2.6), we find that 
equations (2.16) and fl2~"2]) reduce to 
dx 

dn 

dy 



dx 

— = (l-x)[(3j-2)x + y], (2.6a) 



- y[2 + y + (37 - 2)x] + bx ri ~ m y + 9ax ri . (2.66) 

Ct/i L 

Also, from the Friedmann equation, (2.1c), we obtain 

l-x = -^ 3 R9- 2 . (2.7) 

Thus, the line x = 1 divides the phase space into three invariant sets, x < 1, x = 1, 
and x > 1. If x = 1, then the model is necessarily a flat FRW model, if x < 1 then the 
model is necessarily an open FRW model, and if x > 1 the model is necessarily a closed 
FRW model. 

The equilibrium points of the above system all represent self-similar cosmological 
models, except in the case 7 = 3£o- If 7 ^ S(q, the behaviour of the equations of state, 
equation (2.6), at the equilibrium points, is independent of the parameters m and n; 
namely the behaviour is 

Cocpa, an d /5 oc p _1 . (2.8) 

Therefore, natural choices for m and T\ are respectively 1/2, 1. We note that in the 
exceptional case 7 = 3£o, there is a singular point {x — l,y — —37} which represents a 
de Sitter solution and is not self-similar. (This is also the case in the Eckart theory as 
was analyzed by Coley and van den Hoogen M.) 

To further motivate the choice of the parameter ri, we consider the velocity of a 



viscous pulse in the fluid |L5 



1 y /2 

T ' (2 ' 9) 

where v = 1 corresponds to the speed of light. Using (|2.4| ) and equations (2.6), we 
obtain 

v = (ax 1 ' 1 - 1 ) 1 / 2 . (2.10) 
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Now, if r% — 1 then not only do we obtain the correct asymptotic behaviour of the 
equation of state for the quantity /3 but we are also allowed to choose a < 1 since 
then the velocity of a viscous pulse is less than the velocity of light for any value of the 
density parameter x. Thus in the remainder of this analysis we shall choose r\ — l. In 
order for the system of differential equations (2.7) to remain continuous everywhere, we 
also assume m < 1. 



3. Qualitative Analysis 

3.1. m — r\ — l 

We now study the specific case when m — r\ — 1. In this case there are three singular 
points, 

(0,0), (l,y-), and (l,y+), (3.1) 

where 

_ 6-3 7 yg - 3 7 ) 2 + 36a + b - 3 7 , y/(b - 3 7 ) 2 + 36a 
1/ = — g 7, an d V = + ^ • (3-2) 

The point (0, 0) has eigenvalues 
3 7 -4 + 6 1 r- — — - — — 3 7 -4 + 6 1 



-v/(6- 3 7 ) 2 + 36a, -L- + - y/(b - 3 7 ) 2 + 36a. (3.3) 

This point is either a saddle or a source depending on the value of the parameter, B\\ 

B l = (2-6)(37-2) + 9a. (3.4) 

If 5i > 0, then the point is a saddle point and if B\ < 0, then the point is a source. If 
Bi = (the bifurcation value), then the point is degenerate (discussed later). 
The point (l,y~) has eigenvalues 

y/(b - 3 7 ) 2 + 36a, - 37 ~ 2 4 + & + l - V(6 - 3 7 ) 2 + 36a. (3.5) 

If 5i < 0, then the point (l,y~) is a saddle point and if Bi > 0, then the point (l,y~) 
is a source. If Si = (the bifurcation value), then the point (l,y~) is degenerate 
(discussed later). 

The singular point (1, y + ) has eigenvalues 

-^(6-37)2 + 36a, - 37 ~ 4 + & - ^ v 7 ^ ~ 3 7 ) 2 + 36a. (3.6) 

This singular point is a sink for 7 > 2/3 (See also Table [1| for details). 

In addition to the invariant set x = 1, there exist two other invariant sets. These 
are straight lines, y = m±x, where 

(b - 3 7 ) ± y/(b - 3 7 ) 2 + 36a 
m± = )! — . (3.7) 
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The invariant line y = m + x passes through the singular points (0,0) and (l,y + ) while 
the line y = m_x passes through the singular points (0, 0), and (1, y~). These invariant 
sets represent the eigendirections at each of the singular points [see also |Appendix A |. 



In order to sketch a complete phase portrait, we also need to calculate the vertical 
isoclines which occur whenever dx/dVL = 0. From (2.7) we can see that this occurs 
either when x = 1 or when y = —(37 — 2)x. This straight line passes through the origin, 
and through the singular point (l,y~) if B\ = 0. If B\ > then the vertical isocline 
has a negative slope which is greater than the slope of the slope of the invariant line 
y = rri-X, [i.e., m_ < —(37 — 2)], and when B\ < the vertical isocline has a negative 
slope which is less than the slope of the invariant line y = m_x [i.e., m_ > —(37 — 2)]. 

To complete the analysis of this model we need to analyze the points at infinity. 
We do this by first converting to polar coordinates and then compactifying the radial 
coordinate. We change to polar coordinates via 

r 2 = x 2 + y 2 and 9 = tan~ x -, (3.8) 

x 

and we derive evolution equations for r and 9. We essentially compactify the phase 
space by changing our radial coordinate r and our time Q as follows, 

and — = (l-f), (3.9) 



1 + r dr 

that is, the plane M 2 is mapped to the interior of the unit circle, with the boundary of 
this circle representing points at infinity of M 2 . We have (for r\ = 1 and general m) 

dv I 

_ = (1 - f)< f(l - f)[(37 - 2) cos 2 9-2 sin 2 9 + (1 + 9a) cos#sin 9] 
dr [ 

-f 2 [(37-2)cos# + sin#] + f 2 ~ m (l - f) m [6sin 2 tfcos 1 -™ j, (3.10) 
^ = (1 - f) 1 9a cos 2 9 - sin 2 9 - 3 7 cos 9 sin 9 + of ^"(l - f) m ~ l sin 9 cos 2 ~ m J. (3.11) 

We easily conclude that if m = 1 (or any m > 0), then the entire circle, f = 1, is singular. 
Therefore, we have a non-isolated set of singular points at infinity. To determine their 
stability we look at the sign of df/dr as f — > 1. In this case we see 



dr 
dr 



(37 - 2) cos 9 + sin 9 (3.12) 



which implies that points above the line y = — (37 — 2)x are repellors, while those points 
which lie below the line are attractors. 

For completeness, we would also like to determine the qualitative behaviour of the 
system at the bifurcation value B\ = where the singular points are (1, 6 — 2) and the 
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line of singular points y = — (3 7 — 2)x. (Note that since Bi = 0, 6 — 2 > 0.) Fortunately 
we are able to completely integrate the equations in this case to find 



16-2 



y\ 



k\l 



x\ 



(3.13) 



where k is an integration constant. We see that all trajectories are straight lines that 
pass through the point (1, 6 — 2). It is straightforward to see that the line of singular 
points are repellors while the point (1, 6 — 2) is an attractor. We are now able to sketch 
complete phase portraits (See Figures |HT| , [B2| and [53]) . 



3.2. m—1/2 and r x — 1 

This is a case of particular interest since it represents the asymptotic behaviour of the 
FRW models for any m and ri (since at the singular points the viscosity coefficient 
behaves like £ oc p 1 / 2 and the relaxation time like /3q oc p _1 ). Note that the physical 
phase space is defined for x > 0, but the system is not differentiable at x — 0. In this 
case there are four singular points, 



(0,0), 



where 



x 



(9a + 2(3 7 -2)) : 



[x,y), (1,2/ ), and (l,y + ), 



y = -(37 - 2)x, 



(3.14) 



(3.15) 



6 2 (3 7 -2) 2 ' 

and y + and y~ are given be equation (|3.2|). 

The dynamical system, (2.7) is not differentiable at the singular point (0,0). We 
can circumvent this problem by changing variables to u 2 = x and a new time variable r 
defined by dVt/dr = u. The system then becomes 

^ = (l- M 2 )[(3 7 -2) M 2 + 2/ ], (3.16a) 

j- = u[9au 2 -2y-y 2 - (37 - 2)u 2 + byu\. (3.166) 

In terms of the new variables the system is differentiable at the point u = 0, y = 0, but 
one of the eigenvalues is zero and hence the point is not hyperbolic. Therefore in order 
to determine the stability of the point we change to polar coordinates, and find that the 
point has some saddle-like properties; however, the true determination of the stability is 
difficult. [We investigate the nature of this singular point numerically — see [Appendix] 



The singular point (x, y) has eigenvalues 
(3 7 - 2) 2 + 9a 



2(3 7 - 2) 




-2(1 -x)[9a + 2(37-2)]. 



(3.17) 
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This singular point varies both its position in phase space with its stability depending 
upon whether x is less than, equal to or greater than one. If B\ > 0, then x > 1 and 
the point is a saddle point. If B\ < 0, then x < 1 and the point is a source. Finally, if 
Bi = (the bifurcation value), then x = 1 and the point is degenerate (discussed later). 

The stability of the points (l,y~) and (l,y + ) is the same as in the privious case, 
see equations ( |3.5| ) and ( |3.6| ) for their eigenvalues and the corresponding text. (See also 
Table [IJ for details). 

The vertical isoclines occur at x = 1 and y = — (37 — 2)x. This straight line is 
easily seen to pass through the origin and the point (x, —(37 — 2)x). If B\ > 0, then 
the vertical isocline lies below the point (1, y~) and if B x < 0, the vertical isocline lies 
above the point (l,y~). Finally if B\ = (the bifurcation value), the vertical isocline 
passes through the point (l,y~). 

From an analysis similar to that in the previous subsection, we conclude that there is 
a non-isolated set of singular points at infinity. Their qualitative behaviour is the same in 
this case as in the previous case; namely, points which lie above the line y = —(37 — 2)x 
are repellors, while those points which lie below the line are attractors. 

At the bifurcation value £>i = 0, the points (1, y~) and (x, — (37— 2)x) come together; 
consequently these points undergo a saddle-node bifurcation as B\ passes through the 
value 0. The singular point is no longer hyperbolic, but the qualitative behaviour near 
the singular point can be determined from the fact that we know the nature of the 
bifurcation. Hence the singular point is a repelling node in one sector and a saddle in 
the others. A complete phase portrait is sketched in Figures p|, IBS] and IBB. 



Table 1. Qualitative nature of the singular points of the dynamical system, (2.7), for 
different values of the parameter Bi a and r\ = 1 (with respect to fi-time). 







(0,0) 


(l,2T) b 


(i,y + ) c 


(x,y) d 


TO = 


1 and Bi < 


source 


saddle 


sink 




TO = 


1 and Bi = 


source e 


source e 


sink 




TO = 


1 and Bi > 


saddle 


source 


sink 




TO = 


1/2 and B x < 


saddle 


saddle 


sink 


source 


TO = 


1/2 B x = 


saddle 


saddle-node 


sink 


saddle-node f 


TO 


1/2 and B x > 


saddle 


source 


sink 


saddle 



a Bx = (6-6)(7-2)+3a. 

b y~ = (b - 3 7 - yj(b - 3 7 ) 2 + 36a) /2. 

c y+ = (b - 3 7 + y/(b - 3 7 ) 2 + 36a) /2. 

d x = (9a + 2(3 7 - 2)) 2 /6 2 (3 7 - 2) 2 , y = -(3 7 - 2)x, 

e These points are part of the non- isolated line singularity y — —(3j — 2)x. 

1 This is the situation when the points (l,y~) and (x,y) coalesce. 
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4. Discussion 

4-1. Exact Solutions 

The exact solution of the Einstein field equations at each of the singular points represent 
the asymptotic solutions (both past and future) of FRW models with a causal viscous 
fluid source. The solution at each of the singular points represents a self similar 
cosmological model except in one isolated case [see the singular point (l,y~)]. 
At the singular point (0, 0) we have (after a re-coordinatization) 

0(t) = St' 1 , a(t) = a a t, 

pit) = o, n(t) = 0, 

which represents the standard vaccuum Milne model. 

The singular point (l,y + ) represents a flat FRW model with a solution (after a 
re-coordinatization) 

6{t) = 3A + t~ 1 , a(t)=a Q t A+ , 

P (t) = 3(A+) 2 r 2 , n(t) = y + {A+)H-\ 

where A + = 2/(37 + V + ) > 0. 

The singular point (l,y~) represents a flat FRW model. If 7 ^ 3( Q then the solution 

is 

9(t) = 3A~(t - Q- 1 , a(t) = a \t - t Q \ A ~ , 

p(t) = 3(A-)\t - t )~\ U(t) = y + {A-f{t - t o y\ 

where A~ = 2/(37 + V~)- (Note that in this case we cannot simply change coordinates 
to remove the constants of integration.) The sign of A~ depends on the sign of 7 — 3£ . 
If 7 > 3( then A~ > 0, and if 7 < 3£ then A~ < 0. Thus if A~ < then 9 is positive 
only in the interval < t < t a and hence we can see that after a finite time t a 9, p, and a 
all approach infinity. (We will see later that the WEC is violated in this case). If A~ > 
then we can re-coordinatize the time t so as to remove the constant of integration, t , 
and the absolute value signs in the solution for a(t). If 7 = 3£ Q then A~ = and the 
solution is the de Sitter model with (after a re-coordinatization) 

9{t) = 3H , a{t) = a e H °\ 

pit) = 3H \ Ii[t) = H 2 y~. 

This exceptional solution is the only one that is not self-similar. It can be noted here 
that this is precisely the same situation that occurred in the Eckart models studied in 
i and Q. 

The singular point (x, y) represents either an open, flat, or closed model depending 
on the value of the parameter B\. The solution in all cases is (after a re-coordinatization) 
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Q(t) = St' 1 , a(t) = a Q t 

pit) = 3xt~ 2 , n(t) = yr 



4-2. Energy Conditions 

The weak energy condition (WEC) states that T ab W a W b > for any timelike vector W a 



2j| . In the model under investigation this inequality reduces to p > and p+p + U > 0. 



Assuming that p = (7 — l)p from here on, the WEC in dimensionless variables becomes 



x > and y > —3^/x. (4-1) 

The dominant energy condition (DEC) states that for every timelike W a , 
T a bW a W b > and T ab W a is non-spacelike [M. Here this inequality reduces to p > 
and — p < p + II < p which when transformed to dimensionless variables becomes 

WEC and y < 3(7 - 2)x. (4.2) 

The strong energy condition (SEC) states that T ab W a W b - \T a a W b W b > 0. 
Here this inequality reduces to p > and p + 3p + 3n > which when transformed to 
dimensionless variables becomes 

WEC and y > -(37 - 2)x. (4.3) 

If we assume that the WEC is satisfied throughout the evolution of these models 
then we find that there are five distinct situations. If 7 < 3( Q , then B\ > and the line 
y = — 372 intersects the line x — 1 at a point y > y~ . If 7 = 3£ , then B\ > and the 
line y = S'jx intersects the line x — 1 at the point y = y~ . If 7 > 3( Q then Bi can be 
of any sign or zero, but the line y = — 3 / ~yx intersects the line x = 1 at a point y < y~. 
If the WEC condition is satisfied throughout the evolution of these models then the 
possible asymptotic behaviour of the models is greatly restricted. 

4-3. Asymptotic Behaviour 

The qualitative behaviour depends on the values of Bi and m. If the parameter m is 
different from unity then there is an additional singular point. This property is also 
present in the Eckart models studied by Burd and Coley Q and Coley and van den 
Hoogen [|l|. The value m = 1 corresponds to the case when the dynamical system, 
(2.7), is polynomial (the only other value of m that exhibits this property is m — 0). 
The value m — 1/2 is of particular interest as it represents the asymptotic behaviour 
of all the viscous fluid FRW models, and also, this is the case when the equation of 
state for £ is independent of 9 (i.e., ( oc p 1 / 2 ). The parameter, Bi, plays a role similar 
to the parameter 9( a — (37 — 2) found in both Burd and Coley || and Coley and van 
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den Hoogen JT|. The value of the parameter Bi determines the stability and global 
behaviour of the system. 

One of the goals of this paper is to determine the generic behaviour of the system 
of equations (2.7). Using the above energy conditions, and in particular the WEC, 
and the phase-portraits (Figures |B 1|-|B6|) , we can determine the generic and exceptional 



behaviour of all the viscous fluid models satisfying the WEC. We are primarily interested 
in the generic asymptotic behaviour of the FRW model with viscosity: If we consider 
the dynamical system, (2.7), as X = F(X) where X = (x,y,a,b,j) C 1R 5 , x, y are the 
variables, and a, b, 7 are the free parameters then generic behaviour occurs in sets of 
non-zero measure with respect to the set M. 5 (except for the flat models in which case 
the state space is a subset of 1R 4 ). For example, the case B\ = is a set of measure zero 
with respect to the set M 5 . All behaviour is summarized in Table 0. 

Typically, if B\ < 0, then the open models evolve from the big-bang visco-elastic 
singularity at (1, y + ) and evolve to the Milne model at (0, 0) [if m — 1] or the non-vacuum 
open model at (x,y) [if m = 1/2]. If B\ < 0, then the closed models evolve from the 
big-bang visco-elastic singularity at (l,y + ) to points at infinity. These particular points 
at infinity correspond to the points where 6 = (point of maximum expansion) and the 
various dimensionless variables breakdown. 

Typical behaviour of models with B x > depends upon the sign of 7— 3( Q . If 7 < 3( a , 
then all trajectories for the open models will violate the WEC and if 7 > 3( Q , then the 
open models evolve from the big-bang visco-elastic singularity at (l,y + ), become open 
models and then evolve towards an inflationary flat FRW model at the point (1, y~). A 
visco-elastic singularity is a singularity in which a signifigant portion of the initial total 
energy density is viscous elastic energy, that is II ^> 1. Concerning the closed models 
when B\ > 0, if 7 < 3( Q then models evolve from the big-bang visco-elastic singularity 
at (l,y + ) to points at infinity. However, if 7 > 3( OJ then the closed models again evolve 
from the big-bang visco-elastic singularity at (1, y + ) but now have two different typical 
behaviours. There is a class of models which approach points at infinity and do not 
inflate and there is a class of models which evolve towards the inflationary flat FRW 
model at the point (l,y~). 

The flat FRW models consist of a subset of measure zero of the total state space 
M 5 . However, the flat models are of a special interest in that the flat models represent 
the past asymptotic behaviour of both the open and closed models. If B\ < 0, then the 
flat models evolve from the visco-elastic singularity at (1, y + ) to points at infinity or to 
the flat model located at (l,y~). If B\ > 0, and 7 < 3( a , then the flat models evolve 
from the big-bang visco-elastic singularity at (l,y + ) to points at infinity. And if B\ > 
and 7 > 3(o, then models evolve from the visco-elastic singularity at (1, y + ) to points 
at infinity (non-inflationary) or to the inflationary model at the point (l,y—). 

Note, that if the WEC is dropped (i.e., let 7 < 3£ ) then a class of very interesting 



14 



Table 2. Asymptotic behaviour of the FRW models with bulk viscosity satisfying the 
WEC with n = 1. 



parameters m 



models generic behaviour 



exceptional behaviour 3 



7 < 3Co, Bi > m = 1,1/2 open 

fiat 
closed 

7>3Co, Bi > m = 1,1/2 open 

fiat 



7 > 3Co, B x = m=l 



closed 

open 
flat 

closed 



m=l/2 open 
flat 



7 > 3Co, Si < m=l 



m=l/2 



closed 

open 

flat 



closed 

open 

flat 



closed 



(1,2/+ 
0,2/+ 
0,2/+ 

(i,2/+ 
(i,2/+ 
(i,2/+ 
(i,2/+ 
(i,y + 
(i,2/+ 
(i,2/+ 
0,2/+ 
0,2/ + 
0,2/+ 
(i,2/+ 
(i,2/+ 
(i,2/+ 
0,2/+ 
(i,2/+ 
(i,2/+ 
(i,2/+ 
(i,y + 



oo 
oo 

(1,2/-) 

oo 

(1,2/-) 

oo 

(1,2/-) 

(x , -(37 - 2)x ) 
oo 

(1,2/-) 

oo 

o , -(37 - 2)x ) 

(1,2/-) 

oo 

(1,2/-) 

oo 

(0,0) 

oo 

(1,2/-) 

oo 

0,27) 
oo 

(1,2/-) 

oo 



(1,2/+) -(0,0) 

0,27) -> oo, (l,y+) ->■ (x,y) b 
(1,2/+) -(0,0), (0,0) -(l.jf 



0,2/) — (1,2/ ) 

0,27) — oo, (1,2/ + ) — 0,y) b 



(1,2/+) -(0,0) (0,0)->(l,y-) 



(1,2/ ) — 00 
(l,y-) -,(0,0) 



(1,2/ ) -> oo 

(1,2/+) -(0,0), (I,?/") -0,27) 



(1,2/ ) -> oo 



a These are exceptional trajectories and consequently do not represent typical or generic 
behaviour. 

b In the case m = 1/2. 



models occurs. There will exist models that will evolve from the visco-elastic singularity 
at (1, y + ) with 6 > and 6 < start inflating at some point U and at a finite time after U 
will start expanding at increasing rates, that is 9 > 0, and will eventually evolve towards 
the point (1, y~). (This is the special case mentioned in the previous subsection.) What 
this means in terms of the open and flat models is that they will expand with decreasing 
rates of expansion, start to inflate, and then continue to expand with increasing rates 
of expansion. For the closed models, the models will expand with decreasing rates 
of expansion, start to inflate, and then continue to expand with increasing rates of 
expansion, these models will not recollapse. 
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5. Conclusion 



The only models that can possibly satisfy the WEC and inflate are those models with 
7 > 3£ and B\ > 0. Therefore, we can conclude that bulk-viscous inflation is possible 
in the truncated Israel- Stewart models. However, in the models studied by Hiscock and 
Salmonson || inflation did not occur (note, the equations of state assumed in || are 
derived from assuming that the universe could be modelled as a Boltzmann gas), while 
inflation does occur in the models studied by Zakari and Jou [pjl who utilized different 
equations of state. Furthermore, Maartens []22]| has also analyzed models arising from 
the Israel-Stewart-Hiscock theory of irreversible thermodynamics. Maartens assumes 
an equation of state for the temperature T of the from T = (3e~ rHot and finds that the 
inflationary attractor is unstable in the case m > 1/2 |p2| . In our truncated model we 
choose dimensionless equations of state and find that inflation is sometimes possible. 
The question of which equations of state are most appropriate remains unanswered, and 
clearly the possibility of inflation depends critically upon the equations of state utilized 

E3- 

This work improves over previous work on viscous cosmology using the non-causal 
and unstable first order thermodynamics of Eckart || and differs from the work of 



Belinskii et al. |T2| in that dimensionless equations of state are utilized. From the 



previous discussion we can conclude that the visco-elastic singularity at the point (1, y + ) 
is a dominant feature in our truncated models. This singular point remains the typical 
past asymptotic attractor for various values of the parameters m and B\. This agrees 
with the results of Belinskii et al. fll2 |. The future asymptotic behaviour depends upon 
both the values of m and B\. If B\ < 0, then the open models tend to the Milne model 
at (0, 0) [m = 1] or to the open model at (x, y) [m = 1/2], and if B\ > 0, then the open 
models tend to the inflationary model at the point (l,y~) or are unphysical. If B\ < 0, 
then the closed models tend to points at infinity, and if B x > 0, then the closed models 
tend to the inflationary model at the point (l,y~). The future asymptotic behaviour of 
the flat models is that they either tend to points at infinity or tend to the point (l,y~), 



in agreement with the exact solution given in [14 



Belinskii et al. [12] utilized the physical variables 9, p, and II in their analysis and 



assumed non-dimensionless equations of state, and they found a singularity in which the 
expansion was zero but the metric coefficients were neither infinite nor zero; the authors 
passed over this observation stating that in a more realistic theory this undesirable 
asymptotic behaviour would not occur. We note that this behaviour did not occur in 
our analysis in which dimensionless variables and and a set of dimensionless equations 
of state were utilized. 

The behaviour of the Eckart models in Burd and Coley [|J with 9£o — (37 — 2) > 
is very similar to the behaviour of the truncated Israel-Stewart models studied here in 
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the case B 1 < 0. This result also agrees with the conclusions of Zakari and Jou |15 
However, when B x > various new possibilities can occur; for instance, there exist 
open and closed models that asymptotically approach a flat FRW model both to the 
past and to the future. Interestingly enough this is also the case in which the future 
asymptotic endpoint is an inflationary attractor. This type of behaviour does not occur 
in the Eckart theory. 

While this work (and the work of Belinskii et al. |T2|) employs a causal and stable 
second order relativistic theory of thermodynamics, only the 'truncated' version of the 
theory has been utilized p2| , rather than the full theory of Israel-Stewart-Hiscock 



P, |T|, [LI], |23| . However, the 'truncated' equations can give rise to very different 
behaviour than the full equations; indeed, Maartens argues that generally the 
'truncated' theory will lead to pathological behaviour, particularly in the behaviour 
of the temperature. It is possible to reconcile the truncated and full theories, at 
least near to equilibrium, by using a "generalized" non-equilibrium temperature and 



pressure defined via a modified Gibbs equation in which dissipation enters explicitly |24 



Moreover, it is expected that the truncated theory studied here is applicable at least in 
the early universe. We have neglected the divergence terms here as a first approximation. 
If one includes the divergence terms then the system of equations describing the model 
requires an additional equation of state for the temperature T. A reasonable assumption 
in this case might be to assume a dimensionless equation of state for T. 

Notwithstanding these comments it is clear that the next step in the study of 
dissipative processes in the universe is to utilize the full (non-truncated) causal theory 
|8]. [□], |T|, ^3|, and this will be done in subsequent work. In particular, we shall find that 
the work here will be useful in understanding the full theory and is thus a necessary 
first step in the analysis. In addition, anisotropic Bianchi type V models which include 
shear viscosity and heat conduction will also be investigated. 
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Appendix A. First Integrals 



We will use Darboux's theorem |25] to find an algebraic first integral of the system in 
the case m = r x = 1 by first finding a number of a algebraic invariant manifolds. An 
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algebraic invariant manifold, Qi = 0, is a manifold such that Qi = TiQi, where is a 
polynomial. The following are invariant manifolds of the system: 

Qx — x-1, 
Q 2 = y - m_x } 
Q 3 = y - m+ x } 

where 



m± = ( t -3 7 )± > /(6-^ + 3ga - (A1) 
Calculating Qi = r^Qi, we find 

n = -[y + (37 - 2)x], 

r 2 = ~[y + (37 - 2)x + (2 - 6 + m_)], 
^3 = - [2/ + (37 - 2)z + (2 - 6 + m+)]. 

Using Darboux's Theorem, an algebraic first integral Q can be found by setting 
Q — Q^Q^Q^ an d then determining what values of ctj satisfy the equation Q — 0. 
Solving the resulting algebraic system we find the following algebraic first integral of 
the dynamical system, (2.7), in the case m = r± = 1: 

Q = (x-l) Ql (y-m_x) a2 (?/-m + x) a3 = AT (constant) (A.2) 

where a,\ is a free parameter and a.2 and «3 must satisfy 

1 6 + 37-4 \ 

" 2 1 ~2 ~ ^(fe - 37)2 + 36^ J Ql ' 

1 6 + 37-4 \ 

Q3= '"2 + 7^fTSj ai - 

This first integral determines the integral curves of the phase portraits in Figures 
Bl , [B2] and [B3| , where the value determines which integral curve(s) is being described. 
For example, if K = 0, the integral curves are x = 1 and y = m±x. Also, we can see 
that if 6 = 4 — 37, a,\ = —2 and K — 1 then the integral curve describes an ellipse; 
however, these closed curves necessarily pass through the points (l,y + ) and (l,y~), 
thereby nullifying the possible existence of closed orbits. 



Appendix B. Numerical Analysis 

If m = 1/2, the dynamical system (2.7) is only defined for x > 0, but more importantly 
it is not differentiable at x — 0. In this section we will use numerical techniques to 
analyze the integral curves in the neighborhood of the singular point (0, 0) in the case 
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m = 1/2. The integration and plotting was done using Maple V release 3. From the 
qualitative analysis we find that the behaviour depends on the parameter B\. In the 
first of these two plots we choose 7 = 1, a = 1/9, and b = 4, so that B\ = — 1 < (see 
Figure P7|) . In the second plot we choose 7 = 1, a — 1/9, and b = 2 so that B\ — 1 > 
(see Figure BS). From the numerical plots we can conclude that the point (0,0) has 
a saddle-point like nature, in agreement with preliminary remarks made in the text in 
section 
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Figure captions 



Figure Bl. The phase portrait describes the qualitative behavior of the FRW models 
with bulk viscous pressure in the case m = r\ = 1 and B\ < 0. The arrows in the 
figure denote increasing fl-time (fl — > oo) or decreasing t-time (t — > + ). 



Figure B2. The phase portrait describes the qualitative behavior of the FRW models 
with bulk viscous pressure in the case m = n = 1 and B\ = 0. The arrows in the 
figure denote increasing fl-time (fl — > oo) or decreasing t-time (t — > + ). 



Figure B3. The phase portrait describes the qualitative behavior of the FRW models 
with bulk viscous pressure in the case m = r\ = 1 and B\ > 0. The arrows in the 
figure denote increasing fi-time (fl — > oo) or decreasing t-time (t — > + ). 



Figure B4. The phase portrait describes the qualitative behavior of the FRW models 
with bulk viscous pressure in the case m = 1/2 and r\ = 1 with _Bi < 0. The arrows 
in the figure denote increasing fl-time (fl — > oo) or decreasing i-time (f — > + ). 
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Figure captions 



Figure B5. The phase portrait describes the qualitative behavior of the FRW models 
with bulk viscous pressure in the case m — 1/2 and r\ = 1 with B\ = 0. The arrows 
in the figure denote increasing fi-time (17 — > oo) or decreasing i-timc (f — > + ). 



Figure B6. The phase portrait describes the qualitative behavior of the FRW models 
with bulk viscous pressure in the case m = 1/2 and n = 1 with i?i > 0. The arrows 
in the figure denote increasing fi-time (Q — ► oo) or decreasing f-time (t — > + ). 



Figure B7. The phase portrait describes the qualitative behavior of the FRW models 
in the case m = 1/2 and B\ = — 1 < 0. 



Figure B8. The phase portrait describes the qualitative behavior of the FRW models 
in the case m = 1/2 and B\ = 1 > 0. 
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